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Abstract 

In many growth processes particles are highly mobile in an active layer at the surface, 
but are relatively immobile once incorporated in the bulk. We study models in which 
atoms are allowed to interact, equilibrate, and order on the surface, but are frozen 
in the bulk. Order parameter correlations in the resulting bulk material are highly 
anisotropic, reflecting its growth history. In a flat (layer by layer) growth mode, 
correlations perpendicular to the growth direction are similar to a two dimensional 
system in equilibrium, while parallel correlations reflect the dynamics of such a 
system. When the growing film is rough, various couplings between height and order 
parameter fluctuations are possible. Such couplings modify the dynamic scaling 
properties of surface roughness, and may also change the critical behavior of the 
order parameter. Even the deterministic growth of the surface profile can result in 
interesting textures for the order parameter. 



1 Introduction 



For many technological applications, high quality films are grown by the pro- 
cess of vapor deposition. The properties of such films can be quite different 
from the same material produced in bulk equilibrium [1,2], reflecting their 
preparation history. For example, during the growth of some binary alloys, 
the deposited atoms are highly mobile on the surface, but relatively immo- 
bile in the bulk [2]. Consequently, the surface fluctuations occurring during 
the growth process are frozen into the bulk. A characteristic feature of such 
(metastable) phases is anisotropic correlations related to the growth direction 
which are absent in bulk equilibrium. 

A number of models for composite film growth have been introduced in the 
past [3-9]. Generally in these models, the probability that an incoming atom 
sticks to a given surface site depends on the state of neighboring sites in the 
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layer below. Once a site is occupied, its state does not change any more, and 
thus the surface configuration becomes frozen in the bulk. Such growth rules 
are equivalent to (stochastic) cellular automata, where each site is updated in 
parallel as a function of the states of its neighbors. Subsequent states of the 
cellular automaton correspond to successive layers in the crystal. 

It is in general not possible to calculate exact correlation functions for such 
(non-equilibrium) growth processes. The exception occurs in special cases 
where the growth rules satisfy a detailed balance condition, relating their sta- 
tionary behavior to an equilibrium system of one lower dimension [10]. How- 
ever, it can be shown that if <i-dimensional probabilistic cellular automata 
with two states, and up-down symmetry, undergo a symmetry breaking, their 
critical behavior is identical to the corresponding Ising model in equilibrium 
[10]. Correlations in time are then equivalent to those generated by Glauber 
dynamics of the Ising system, (d + l)-dimensional crystals grown according 
to the rules of these cellular automata therefore have an order-disorder phase 
transition with correlations perpendicular to the growth direction character- 
ized by the critical exponent is, and those parallel to the growth direction by 
the exponent vz of the d-dimensional Ising model (z being the appropriate 
dynamical critical exponent). 

In Sec. 2, I introduce a model for layer by layer growth of binary films. The 
atoms on the top layer are assumed to equilibrate completely (by surface diffu- 
sion or desorption-resorption mechanisms) before another layer is added [11]. 
Such an assumption is realistic only if the growth rate is much slower than 
characteristic equilibration times of the surface layer. The model satisfies de- 
tailed balance, and can therefore be analyzed with methods from equilibrium 
statistical physics. This discrete model is then used to justify a continuum 
formulation to the problem which is identical to the time dependent Landau- 
Ginzburg equation for model A dynamics [12]. 

In general, a layer by layer growth mode is unstable, and the role of surface 
roughness is explored in Sec. 3. The Kardar, Parisi, Zhang (KPZ) equation 
describes the dynamic fluctuations in the height of an amorphous surface [13]. 
The interplay between roughness and ordering phenomena is then considered 
by introducing simple equations that couple fluctuations in height and the 
order parameter. Long range correlations occur at the critical point for the 
onset of ordering in the surface binary mixture. This in turn leads to greater 
roughness fluctuations, whose scaling can be explored perturbatively around 
d — 4 dimensions [14]. 

While the order parameter for binary deposition is a scalar, we can more 
generally examine the case of a continuous order parameter formed on the 
surface layer. There are soft modes associated with such continuous symme- 
try breaking whose coupling to height fluctuations are explored in Sec. 4. In 
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particular, the deterministic relaxation of the order parameter on an initially 
rough surface can in fact be described exactly through a generalized Cole-Hopf 
transformation. Interestingly, the relaxation process is super-diffusive and oc- 
curs through coarsening of domains (separated by sharp domain walls) on 
surface mounds. 



2 Layer by Layer Growth 



2.1 Discrete model 



Binary growth is modelled by two kinds of atoms, A and B, which occupy the 
sites of a d+ 1 dimensional hypercubic lattice. Let €aa, £ab, and €bb denote 
the interaction energies between neighboring atoms of types AA, AB, and BB 
respectively. When each layer has ./V sites, there are 2^ possible configurations 
for a layer. The energy cost for adding a new layer of configuration 7 on top of 
one in configuration a is the sum of the internal energy E 1 of the new layer, 
and the interaction energy V ai with the previous layer. These energies are the 
sums of all local bonds between nearest neighbors ij within the new layer, 
and between the two layers, respectively. In addition, E 1 contains a chemical 
potential HaNa + HbNb related to the partial pressures of A and B atoms in 
the gas phase. 

Assuming that the top layer is in thermal equilibrium, the conditional proba- 
bility that it is in configuration 7, given configuration a for the layer below, 
is 

exp [— gQKy + Very)] m 

7 ° E S e X p[-P{E 6 + V aB )y [) 



where T = (ksP)^ 1 is the temperature at which the crystal is grown. After 
adding many layers, the steady-state probability for finding a configuration 7 
is determined by the stationarity condition 

p 7 = ]TwvP a , (2) 

a 



which has the solution 



7 ex P [ — 

a ~ E^exp [-/3(E S + E„ + V Sv )} = E^exp [-(3H 5v \ ' U 



The above expression is the equilibrium probability for the top layer of a 
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two-layer system, obtained after summing over the states of the bottom layer. 
Transverse correlation functions (i.e. perpendicular to the growth direction) 
are therefore exactly the same as correlation functions in a two-layer system. 

From Eqs. (1) and (3) it follows that the system satisfies detailed balance, i.e., 



Thus, beyond a transient thickness, the crystal looks the same along or against 
the growth direction, and the sequence of layers corresponds to time evolution 
of thermodynamic equilibrium states. This generalizes previous results for cel- 
lular automata, which are obtained by setting the in-plane interactions E a to 
zero. As in such cellular automata, the (d+ l)-dimensional system has trans- 
verse properties like rf-dimensional models. In particular, phase transitions 
occur at the same temperature as for a rf-dimensional two-layer system. 

Generalizing the model, by allowing several layers at the surface to equilibrate, 
is straightforward. To mimic the large energy of the impinging particles, as well 
as their modified environment, we can assign each of the top £ layers from the 
surface a different temperature, through scaled interaction energies depending 
on its depth. The probability that a layer with configuration 7 follows one in 
configuration a in the bulk is obtained by considering the layer at the moment 
when it is the £th layer from the top, i.e. immediately before its configuration is 
frozen. Denoting the configuration of the first £—1 layers by C 7 and their energy 
(including the coupling to the £th layer, and different interaction constants in 
the different layers) by E(C 7 ), the conditional probabilities W ia can be written 
as 



Following the approach for the case £ = 1, we can show that the set of weights 



describe a stationary state. It is easy to verify that this stationary solution 
satisfies detailed balance. The stationary state corresponds to an equilibrium 
Hamiltonian with 2£ layers, with interactions which depend on the distance 
from the closest surface. The top (or the bottom) layer describes the deposited 
surface, while the middle (£ or £ + 1) layers describe transverse correlations 
in the bulk. While the correlations parallel to the growth direction are more 
complicated, the general conclusions for £ = 1 remain valid. 




W ia P a . 



(4) 




Ec 7 exp{-/?[E 7 + \/ a7 + E(C 7 )]} 
£«A ex P i-Pl E s + v oS + E(Cs)]} ' 



E 7A ,c Q exp {-P[E a + E 7 + V ai + ggj + E(C a )]} 
Es,v,c s ,c„ exp {-/3[E S + E v + V Sv + E(C S ) + E{C„)]} ' 
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2.2 Continuum formulation 



In the above discrete model, we can use an Ising variable <jj = ±1 to indicate 
if site i is occupied by atom A or B. Close to the critical point, density fluctua- 
tions occur over long distances and universal properties are better captured by 
considering a coarse-grained order parameter m(x, t). Here x labels the d di- 
rections transverse to growth, while t which indicates time is also proportional 
to the coordinate parallel to the growth direction. Hence m(x, t) encodes the 
time history of the growth process. From the exact solution of the discrete 
problem, we know that the behavior of these configuration is equivalent to 
the time evolution of a <i-dimensional system at equilibrium. In the contin- 
uum limit, the latter is described by the time-dependent Landau Ginzburg 
equation [12] 

d t m = KV 2 m + rm — urn? + r] m (x, t), (5) 

where ^ m (x, t) is a random noise of zero mean, whose variance is proportional 
to the growth temperature. 

Away from the critical point at r = r c , fluctuations in m decay over a trans- 
verse correlation length £, and a longitudinal correlation 'time' £ 2 . At the 
critical point itself, there is no intrinsic scale, and correlations decay as 

(m(x, t)m(x', t')) = | \r^-g m ( • (6) 

|x — x'| 2Xm \l x — x | / 

In dimensions d > 4, criticality occurs for r = u = (the diffusion equation), 
leading to z = 2 and Xm = (2— d)/2. On approaching criticality, the correlation 
length diverges as £ oc |r — r c \~ u , with v — 1/2. For d < 4, the nonlinear term 
um 3 is relevant, and the exponents can be calculated perturbatively[12] in 
e = 4 - d. 



3 Rough Growth 

3.1 Dynamic roughening 

The layer by layer growth mode cannot be maintained indefinitely, and the 
surface eventually becomes rough[15]. Let us denote the height of the surface 
at location x at time t by a function h(x,t). There is considerable evidence 
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from simulations (and some experiments) that the resulting surfaces exhibit 
self-affine fluctuations [16], well described by the continuum equation[13] 

d t h = v + uV 2 h + ^(Vh) 2 +r} h (x.,t). (7) 

The first three terms in this equation correspond to the average deposition 
flux, relaxation by evaporation, and lateral growth, respectively The last term 
represents the fluctuations in the deposition flux, and has zero average. 

The self-affine fluctuations in the surface height can be described by dynamic 
scaling exponents Xh an d z. For example, the averaged two point correlations 
behave as 

([Mx, t) - M< t')] 2 ) = |x - A 2xh 9K ( • ( 8 ) 

The linear equation for A = gives diffusive exponents Xh — (2 — d)/2 and 
z — 2. Any nonlinearity is relevant in d < 2, while sufficiently large A is 
required in d > 2 to produce a rough phase (xh > 0). 

3.2 Coupling growth and ordering 

There are few studies of the interplay between fluctuations in height and 
the order parameter. Some numerical simulations have incorporated both ele- 
ments: As a model for diamond growth, Capraro and Bar-Yam[17] introduced 
a variant of ballistic deposition which exhibits sublattice ordering. Kotrla and 
Predota[9] have examined binary deposition in 1+1 dimensions, resulting in 
domains with rough surfaces. In a recent work with Barbara Drossel[14], we 
took an analytical approach to this problem. 

The starting point is the continuum Eqs. (5,7) describing the order parameter 
ra(x, t), and height h(x., t), fluctuations. To these equations we added all terms 
consistent with the symmetries of the problem. The lowest order (potentially 
relevant) terms result in the following pair of coupled differential equations 

' d t h= vV 2 h + \(yhf + (h ~%rn 2 
d t m = KV 2 m + rm — um 3 + ( m + aVh ■ Vm +bmV 2 h + |m {Vh) 2 

(Note that these equations satisfy the symmetry m \— > —m.) Fluctuations of 
the surface are modified by coupling to the order parameter, through the term 
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proportional to am 2 . There are also three coupling constants a, b, and c, which 
modify the order parameter fluctuations due to coupling to h. 

As long as the binary mixture is disordered (r > r c ), fluctuations in m and 
hence m 2 , are short-ranged, and am 2 acts as another source of white noise. 
The surface fluctuations should thus scale with the standard KPZ exponents. 
However, the range of correlations increases as r — > r c and £ ~ \r — r c \ ~ u — > oo. 
This modifies (most likely increases) the overall amplitude of surface rough- 
ness, and height fluctuations over a scale L behave as 

y/(Sh*(L,r)) = &-x»Lx»g (L/0 , (10) 
where Xh * s the roughness exponent at criticality, which is discussed next. 

3.3 Critical roughness 

Under a change of scale x i— > fex, t i— > ft, i— > 6 x,i /i, and m i— > b Xm m, the 
non-linear coefficients in Eqs. (9) scales as rr i— > 6 y:E a;, with 

2/a = 2/a = y& = Xft + ^ - 2, y c = 2x^ + 2-2, y a = 2x m - Xh + z. (11) 

The critical point in dimensions d > 4 occurs at r = u = 0. The linear diffusion 
equations at this point result in the bare field dimensions Xh = Xm = (2 — d)/2. 
Taking account of the non-linearities, we observe the following behaviors. 

• d > 6 : All non-linearities are (perturbatively) irrelevant; the surface is 
smooth, and the order parameter goes through a classical phase transition. 

• 4 < d < 6 : The leading non-linearity is the term am 2 describing the 
correlated noise acting on the surface height, with (z = 2) 

yl = A-d-x c h = Z-d/2. (12) 

In these dimensions, the correlated noise is more relevant than the white 
noise from the flux variations [18]. The correct result can in fact be obtained 
simply by setting y a to zero, leading to critical height fluctuations with 

Xl = ^d> 2 -^ . (13) 

Note that while the roughness exponent is larger than its bare value, it 
is still negative. The scaling of the order parameter is not modified, and 
Xm = (2-d)/2. 
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• d < 4 : When the roughness exponent is positive, all the couplings A, a, 
b, and c become relevant. Also in d < 4, the critical point of the Landau- 
Ginzburg model is no longer at r = u = 0, and a full renormalization group 
(RG) study is called for [14]. Ignoring the feedback from height fluctuations 
to the order parameter, we find to leading order an RG equation of the form 



where e = 4 — d, and C is a positive constant. 

There is a fixed point at aX = e/C, with roughness exponent xt = 0- m 
d > 4, this is an unstable fixed point governing a transition between flat 
(xl — 4 — d < 0) and rough phases (occurring for aX < — (d — 4)/C). For 
g? < 4, this fixed point is stable and attracts all points with aX > 0. Negative 
values of aX flow to a rough phase which is not perturbatively accessible. 
Including all non-linearities in the equation for m complicates the analysis, 
but we did not find a fixed point whose critical behavior is different from 
the ordinary Landau-Ginzburg model (at least to lowest order). 



4 Continuous Order 

4-1 Stochastic evolution 

The situation on the ordered side of the phase transition is more complex. 
The analogy to the dynamics of the lower- dimensional system suggests that the 
leading process is the gradual coarsening of the ordered domains. Such domains 
would then appear as cone-shaped columns in the bulk film, a reasonably 
common feature of growth textures. However, more work is necessary to verify 
and quantify this picture. 

Another interesting situation is when the symmetry breaking involves a con- 
tinuous, rather than a discrete (Ising like), order parameter. For example, we 
may consider deposition of spins which can realign on the surface but are 
frozen in the bulk. More interestingly, the growth of crystals involves transla- 
tional and orientational symmetry breakings in the plane. In the simplest case 
of a vector order parameter, we can simply generalize Eqs. (9) by replacing 
the scalar m with an n-component vector m(x, t). While the discussion of 
critical roughening is not significantly modified from the Ising case (n = 1), 
new issues arise pertaining to the ordered phase. 

The most common excitations of the broken symmetry phase are soft (Gold- 
stone) modes, which can in principle couple to the surface roughness. The 
simplest example is provide by the XY model (n — 2), where the direction of 




1 d(aX) 
vXj di 



e-C(aX), 



(14) 
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the vector can be described by an angular field 0(x, t). Including the lowest 
order terms which satisfy rotational symmetry leads to the coupled equations 
of motion 

' d t h = vV 2 h + l(Vh) 2 + ( h -f(W) 2 
< - (15) 

d t 9 = KV 2 9 + ( e +aVh ■ V9 



Interestingly, these are precisely the equations proposed in Refs. [19,20] in the 
contexts of moving flux lines and drifting polymers. In these contexts, the 
above equations have been studied by RG analysis (around d — 2), and by 
numerical simulations (in d — 1). In particular, in d = 1 the KPZ exponents 
(Xh = 1/2 and z = 3/2) are recovered for the surface roughness, while the 
angular fluctuations remove any long-range order. Further analysis is again 
necessary for the case d = 2. Specifically, an important aspect of the field 9 
not present in the earlier studies is its angular nature. It could thus include 
vortices which are topological defects. Such defects typically play an important 
role in equilibrium two dimensional systems, and have been recently considered 
in a number of related non-equilibrium situations [21]. 



4-2 Deterministic textures 



It is well-known that the non-linear KPZ equation can be recast as a linear 
diffusion equation through the Cole-Hopf transformation. This transformation 
can in fact be generalized to describe the coupling of the surface height to a 
vector order parameter. Consider a field of unit spins, |s(x, t)\ = 1, and set 



W(x, t) = exp 



A/i(x,t) 
2^ 



s(x,t). 



(16) 



A diffusive equation of the field W(x,t), as 

d t W = v\7 2 W + ^r]h(^t)W 



(17) 



can be recast into the pair of coupled differential equations 



d t h = uV 2 h + f (VhY + (^j s- V 2 s + ^(x, t) 
d t s = v [V 2 s - (s ■ V 2 s ) s] + XVh ■ Vs 



(18) 



Note that the transverse component of V 2 s contributes to d t s, thus ensur- 
ing that the magnitude of s is not changed in time, while the longitudinal 
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component of this quantity couples to the surface height. 



It can be checked easily that for n — 2, the parametrization s — (cos 8, sin 6), 
reduces Eqs. (18) to Eqs. (15) in the special limit of a\ = 4z/ 2 , K — v, a — A, 
and r]e = 0. It is also possible to construct other Cole-Hopf transformations 
for cases when a\ < 0[20]. 

Starting from any arbitrary initial condition at t — 0, the deterministic limit 
(Vh — 0) of Eqs. (17-18) is easily solved using the diffusion kernel, as 



— * 

W(x, t) = exp 



A/i(x,t) 



2u 
d d x.' 



s(x, t) 



d/2 



exp 



(x-x ; ) A/i(x',0) 



2v 



s(x',0). 



(19) 



The saddle-point evaluation of the above integral (formally exact as v — > 0) 
captures the long time behavior of the solution. The surface profile 



/i(x, t) = min 



/i(x',0) - 



x — X 



2Xt 



(20) 



consists of a set of parabolic mounds centered at locations x' = x (x, t) cor- 
responding to high points of the initial surface. Note that the evolution of the 
surface profile in this limit is independent of s. The evolution of spins on the 
other hand is completely controlled by the surface height, and given by 



s(x,t) = s(x o (x,t),0) , 



(21) 



i.e. each of the surface mounds carries the spin of its initial high point! Such 
behavior is quite different from the diffusive evolution of spins in the absence 
of coupling to the surface profile. For self-amne initial surface profiles, the 
relaxation of the height and spins is now both diffusive. Furthermore, the spin 
textures produced by this process are domains separated by sharp domain 
walls, very different from the soft modes and vortices that characterize diffusive 
relaxation. Similar extensions of the Cole-Hopf transformation to matrix order 
parameters are also possible [20], and could for example describe relaxation of 
crystalline substrates. 
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